Quantum criticality of vanadium chains with strong relativistic spin-orbit interaction 
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We study quantum phase transitions induced by the on-site spin-orbit interaction AL ■ S in a toy 
model of vanadium chains. In the A ^- limit, the decoupled spin and orbital sectors are described by 
a Haldane and an Ising chain, respectively. The gapped ground state is composed of a ferro-orbital 
order and a spin liquid with finite correlation lengths. In the opposite limit, strong spin-orbital 
entanglement results in a simultaneous spin and orbital-moment ordering, which can be viewed as 
an orbital liquid. Using a combination of analytical arguments and density-matrix renormalization 
group calculation, we show that an intermediate phase, where the ferro-orbital state is accompanied 
by a spin Neel order, is bounded on both sides by Ising transition lines. Implications for vanadium 
compounds CaV204 and ZnV204 are also discussed. 



Quasi-one-dimensional Mott insulators with strongly 
coupled spin and orbital degrees of freedom have at- 
tracted considerable attention recently. A well-studied 
case is the Kugel-Khomskii Hamiltonian with an SU(2) 
symmetry in both spin and orbital sectorsii This model 
is believed to describe the essential physics of quasi-lD 
compounds Na2Ti2Sb20 and NaV205.^''^ Extensive nu- 
merical and analytical studies have revealed a rich phase 
diagram.i^ Of particular interest is a SU(4) symmetric 
point of the Hamiltonian where the low-energy physics 
is described by a conformal field theory with a central 
charge c = 3, equivalent to a model of three free bosons. 

In this paper, we investigate an ID spin-orbital sys- 
tem which in many aspects is different from the above 
SU(2)x SU(2) model. The interest is partly motivated by 
recent experimental progresses on several vanadates in- 
cluding spinel ZnV204'5''^ and quasi-lD CaV204-^'^ The 
vanadium chains in these compounds are characterized 
by frustrated magnetic interactions, Ising-like orbital ex- 
changes, and a large relativistic spin-orbit (SO) interac- 
tion. The origin of the first two features can be traced 
to the lattice geometry of these compounds: vanadium 
ions in spinel form a three-dimensional pyrochlore lat- 
tice, and in CaV2 04 they are arranged in weakly cou- 
pled zigzag chains (Fig. [IJ . In both structures the 90° 
angle between vanadium-oxygen bonds in a network of 
edge-sharing VOg octahedra makes direct exchange the 
primary mechanism for inter-site spin-orbital interaction. 
An important consequence is that the spin-orbital Hamil- 
tonian only depends on orbitals through the correspond- 
ing projection operators. For example, direct dda ex- 
change takes place along a [110] bond only when one or 
both of the dxy orbitals are occupied, giving rise to an 
Ising-like orbital interaction . ^^'^^ 

Contrary to the anisotropic orbital exchange, magnetic 
interaction governed by the Heisenberg Hamiltonian pre- 
serves the spin-rotational symmetry. The spin exchange 
constant, however, depends on the underlying orbital oc- 
cupations. Combined with geometrical frustration, the 
orbital-dependent spin exchange renders these vanadates 
essentially quasi- ID spin systems. To see this, we note 
that the local dxy orbitals are always occupied due to dis- 



FIG. 1: Vanadium chains in (a) CaV204 and (b) ZnV204. 
The black and white circles denote vanadium and oxygen ions, 
respectively. The V"^"'' ions are arranged in a zigzag chain of 
edge-sharing VOe octahedra in CaV204. On the other hand, 
V'^^ ions in spinel form a pyrochlore lattice, which can be 
viewed as a cross-linking network of vanadium chains. The 
quasi- ID spin-1 chains are highlighted by solid blue lines. 



torted VOg octahedron in both compounds. As a result, 
the largest spin-spin interaction takes place on bonds par- 
allel to [110] and [110] directions in spinels [Fig. [Ub)], 
whereas for CaV204 the dominant spin exchange occurs 
along the two rails of a zigzag chain [Fig. [TJa)]. Further- 
more, couplings between these spin chains are not only 
weak, but also geometrically frustrated J^iH Since V'^^ 
ions in both compounds have spin S — 1, their magnetic 
properties thus can be understood from the viewpoint of 
weakly coupled spin-1 chains (solid lines in Fig. [2]). 

In both vanadates, one of the two electrons of the V'^"'" 
ion always occupies the low-energy d^y state, the other 
one occupies either dyz or d^x orbitals. We introduce a 
pseudospin-1/2 operator r to describe this doublet or- 
bital degeneracy, where t" are the Pauli matrices. We 
choose a basis such that r^ = ±1 corresponds to states 
\yz) and i\zx), respectively. The dominant pseudospin 
interaction is governed by an antiferromagnetic Ising- 
like Hamiltonian on nearest-neighbor bonds connecting 
different spin-1 chainsiiii^ (dashed lines in Fig. [5]); the 
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FIG. 2: Schematic diagram of antiferro-orbital orders in (a) 
CaV204 and (b) TjdSfiOs,- In both cases, the d^y orbital is 
occupied at all sites. The interaction between the remaining 
dzx and d^z orbitals (indicated by red and green symbols, 
respectively) is governed by an antiferromagnetic Ising-like 
interaction on the dashed bonds. Also note the ferro-orbital 
order along the quasi- ID spin-1 chains (solid blue lines). 



relevant orbitals on these bonds are dyz and dzx- 



Due 

to the static and three-dimensional nature of the orbital 
Ising Haniiltonian, the system tends to first develop a 
long-range orbital order upon lowering the temperature. 
Fig. [2] shows a schematic diagram of the 3D antiferro- 
orbital order for the two compounds. It is important to 
note that orbitals on individual spin-1 chains (solid lines 
in Fig. [2]) are ferromagnetically ordered. 

The antiferro-orbital order shown in Fig. ^b), how- 
ever, is incompatible with the observed crystal symmetry 
lAi/amd of tetragonal ZnV204i^ii The discrepancy can 
be attributed to a large SO interaction AL • S of vanadium 
ions. Indeed, the SO term is minimized by a state with 
simultaneous Neel ordering of spins and orbital angular 
momenta J^ii^ The corresponding orbital order consisting 
of complex dyz ± idzx orbitals preserves both the mirror 
inversion m and diamond glide d, and is consistent with 
the experimental data. 

In the absence of SO interaction, the ground state of 
the spin-1 chain is a nondegenerate spin singlet. This 
spin liquid phase, also known as the Haldane phase, must 
be separated from the Neel state favored by a large SO 
coupling by quantum phase transitions. Since a nonzero 
L^ requires the electron be in a complex orbital state 
-j= {\zx) zti\xy)), a fully occupied d^y orbital thus re- 
sults in the vanishing of L^ and L^'iiiii^ The nonzero z 
component of the orbital angular momentum is given by 
L^ = — T^ in our representation. To understand the crit- 
ical behavior of vanadium chains due to the LS coupling, 
we consider the following spin-orbital Hamiltonian: 
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This simple model describes two well-studied ID systems, 
i.e. an S* = 1 Haldane chain and a ferromagnetic Ising 
chain (both J, A' > 0), coupled together by an on-site SO 
interaction A L • S. Note that the eigenstates of r^, \yz) ± 
i\zx), carry an angular momentum L^ = =pl, respectively. 



The Hamiltonian Eq. ([T]) has a U(l)x Z2XZ2 symmetry: 
the spin SU(2) symmetry is reduced to U(l) x Z2 by the 
SO term, whereas an additional Z2 symmetry comes from 
the orbital Ising Hamiltonian. 

If should be noted that the ferro-orbital order along the 
quasi-lD chains is stabilized by the inter-chain antiferro- 
orbital coupling in real compounds. The ferromagnetic 
exchange —K in Eq. ^ thus should be regarded as an 
effective coupling in the mean-field sense. 

Despite its simplicity, the model contains rather rich 
physics. It is easy to see that the first-order correction 
vanishes identically in the ground state of decoupled Hal- 
dane and Ising chains. To have a glimpse of the effects 
of the SO interaction, one needs to go to higher orders 
and examine the elementary excitations of model ([T} . We 
start with the kink excitations of the Ising chain. Kinks, 
or domain walls, are topological defects separating the 
two degenerate ground states of perfectly aligned pse- 
dospins. For classical Ising chains, kinks are static quasi- 
particles with a constant energy 2K. To obtain the quasi- 
particle operators, we first fermionize the Ising chain us- 
ing Jordan- Wigner transformation^^ 

< - n (24c™ - 1) (c„ + 4) , < = 1 - 24c„. (2) 

The kink operator 7^ is obtained after subsequent Fourier 
and Bogoliubov transformations 7^ = UqCq 
where Uq = cos(g/2) and Vq = sm{q/2)<^ 

For spin-1 chain, the lowest excitation above the sin- 
glet ground state is a triplet with a dispersion uJk ~ 
[Ag -f ii^(fc — tt)^]^/^ near the energy minimum. Here 
V — 2.56 J is the spinwave velocity and Aq ~ 0.4J 
is the Haldane gap;^^ To model the low-energy physics 
of spin-1 chain, we follow Ref. |l3 and introduce three 
massive magnons a^ and o^ carrying quantum number 
S" = ±1, and 0, respectively. The spin operator is 
S^^l'n + ^Ek '/'fe e*''"/VI, where 



iVqC'_g, 



I 2v 

l-^TT + k 
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(3) 



and the uniform part Z^ is quadratic in transverse 
magnons a^. . In terms of magnons and kinks, Hamil- 
tonian ^ becomes 



H^Y^^kal^al +2X^7^, 



(4) 
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k,q,q' 



The prime on the summation indicates conservation of 
momentum k — q ± q' . In obtaining the above expres- 
sion, we have neglected the interaction between kinks and 
transverse magnons a^ . 

Assuming X ^ J, K, we employ a perturbation theory 
to examine the properties of quasiparticles in the pres- 
ence of SO coupling. We first evaluate the two-kinks 
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FIG. 3: One-loop corrections to self-energy of (a) magnon 
and (b) kink. The solid and wavy lines denote the kink and 
magnon propagators, respectively. 



bubbles shown in Fig. [3Ja) 



ni(w) 



2\'^K 



uj^-iGK"^' 



(5) 



It is interesting to note that the particle- hole bubble does 
not contribute to Hi. The above expression for ni(a;) di- 
verges as the magnon energy approaches that of a pair 
of kinks, i.e. a; « 4:K. In this regime magnons strongly 
interact with the kinks, and higher-order corrections to 
the interaction vertex have to be taken into account. As- 
suming Lj <^ K and using random-phase approximation 
to compute the magnon self-energy, we obtain a renor- 
malized spin gap 



Ao- 



vX^ 



iAoK' 



(6) 



which decreases with increasing A. At large A, closing 
of the spin gap indicates a phase transition into a spin 
ordered phase characterized by {a-jr) 7^ 0. 

Fig. in^b) shows the one-loop contribution to the self- 
energy of magnons. In the K ^ Aq limit, the second 
term in Fig. [3][b) is negligible compared with the first 
one, we obtain a self-energy 



^0 



(7) 



The energy of the kink excitation given by e, « 2K + 
Si(g) indicates that the kinks become mobile through 
the mediation of virtual magnons. 

The perturbative calculation gives important insight to 
the elementary excitations in the small A regime of the 
spin-orbital model ([1]). In particular, the reduced spin 
and orbital gaps indicate quantum phase transitions at 
finite A. To investigate the nature of the phase transi- 
tions and the properties of possible new phases, we nu- 
merically investigate the spin-orbital model Eq. ^ using 
the infinite-system density-matrix renormalization group 
(DMRG) method.i^ The DMRG calculation is known to 
give a rather accurate description of the ground-state 
properties for ID systems. In our calculation we have 
employed periodic boundary conditions in order to ac- 
commodate the staggered ordering of spins and orbital 
angular momenta. 
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FIG. 4: (a) Order parameters as a function of \/J for a chain 
with K = 0.5J. The dashed lines are fittings to 2D Ising 
transition m ~ |A — Ac|^ , where m is the corresponding 
order parameter, (b) The same plot in log-log scale. 



Noting that L^ — —t^ in our representation, we define 
the following order parameters: 

{<)=r, {Sl)=N{-ir, (LfJ = -C(-ir. (8) 

The numerical results of a spin-orbital chain with ex- 
change constants K = 0.5J are shown in Fig. [D At 
small A, the ground state is characterized by a nonzero 
ferro-orbital order r 7^ 0, while the spin sector is in 
the disordered Haldane phase with M = Q. At crit- 
ical point Aci ~ 0.491J, the linear chain undergoes a 
quantum phase transition into a state with simultaneous 
ordering of staggered spin and orbital-moment, charac- 
terized by nonzero order parameters M and C,, respec- 
tively. The ferro-orbital order r remains finite in the 
intermediate phase. As we further increase A, the sys- 
tem undergoes yet another quantum phase transition at 
Ac2 — 0.657J. This critical point is marked by the melt- 
ing of the ferro-orbital order t. In the A -> 00 limit, both 
order parameters A/", C, approach 1. The ground state of 
the spin-orbital chain consists of alternatively occupied 
states \S^ — ±1) ® \yz) =F i\zx). The orbital occupation 
numbers Uyz = nzx = 1/2 are uniform along the chain. 

Since the order parameters M and r describe respec- 
tively the broken Z2 symmetry of the spin and orbital 
sectors, both critical points Ad and Ac2 are expected to 
be in the 2D Ising universality class. Indeed, by fitting 
the corresponding order parameter m to the Ising scal- 
ing relation m ~ |A — Ad"'^/®, we find agreeable result as 
shown by the dashed lines in Fig. 21 The Ising nature of 
the spin Neel transition can be understood in the limit 
of large orbital gap K ^ J. By integrating out orbitals, 
one obtains an easy-axis spin anisotropy: —DS^, where 
D = \^ /4:K-, the resultant spin gap is reduced in ac- 
cordance with Eq. ^. As demonstrated numerically in 
Ref. [20, the spin-1 chain undergoes an Ising transition 
into a Neel state when D > Dc- 

Above Aci , a nonzero Af exerts an effective (staggered) 
transverse field on the orbitals. By rotating pseudospins 
an angle tt about t" axis on odd-numbered sites, the or- 
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FIG. 5: Phase diagram of spin-orbital model ([T|. The lines 
are guide for the eye. In phase I, the spin sector is in the 
disordered Haldane phase, A/" = C = 0, while the orbitals 
are ferromagnetically ordered r 7^ 0. The ground state of 
phase III is composed of Neel spin order and orbital-moment 
order, which can be viewed as an orbital liquid, i.e. A/" 7^ 
and r = 0. In the intermediate phase II, a spin Neel order 
coexists with the ferro-orbital Ising order. 



bital sector is described by a quantum Ising Hamiltonian: 



-Worbital — J^ 7 , ''"it''"rt-H 



rE<' 



(9) 



where the transverse field F = XN". It is known that the 
Ising chain reaches a critical state at Tc = KJ^ Numeri- 
cally, we obtain an Neel order Af « 0.659 at critical point 
Ac2. The corresponding critical field Tc = \c2J^ ~ 0.52J 
is indeed close to the orbital exchange K — 0.5J used in 
the DMRG calculation. 

Our main results are summarized in the phase diagram 
Fig. \5\ which contains three massive phases separated by 
two Ising transition lines. The spin sector in the J — 
limit is extensively degenerate as each spin could be in 



either \S^ = +1) or \S^ — —1) states, independently; the 
total degeneracy is 2^. After applying a 7r-rotation about 
T^ axis to those pseudospins at sites where S^ = —1, the 
orbital sector is again mapped to a quantum Ising chain 
Eq. ([9]) with the transverse field F = A. Since the orbital 
Ising transition occurs at Tc = K, the phase boundary 
Ac2 hence ends at K — Ac2 on the J = axis. In the small 
J limit, one can estimate Ad for the spin Neel transition 
using the critical condition D — Dc ^ O.OSJ^^ which 
gives JK ex Ad- The region of intermediate phase en- 
closed by boundaries Ad and Ac2 shrinks with increasing 
spin exchange J. At very large J, the two Ising lines 
could merge to form a Gaussian criticality or a first-order 
transition. 

We now discuss implications of our findings to vana- 
dium compounds. As discussed before, the ferro-Ising or- 
der parameter t can also serve as the 3D antiferro-orbital 
order parameter for both vanadates (see Fig. [J). A 
nonzero r thus creates two different orbital chains, hence 
further lowering the crystal symmetry. However, experi- 
ments on both compounds observed a higher symmetry^ 
indicating that vanadium chains in both vanadates are 
likely in the t = orbital liquid state (phase III). Fur- 
thermore, the appearance of finite orbital moment L^ ~ ^ 
antiparallel to spin S"^ at A > Ac2 also explains the re- 
duced vanadium moment /x — (2S^ + L^)^b ~ l/^s ob- 
served experimentallyi^ On a final note, we caution that 
the fermionic description of orbital excitations as kinks 
in an Ising chain is the consequence of using ID approx- 
imation for the orbital system. Furthermore, the Kugel- 
Khomskii-type spin-orbital terms {rfTj) (Si • Sj) intro- 
duce correlations between orbital and spin excitations. A 
detailed description of these 3D spin-orbital excitations 
will be discussed in future publications. 
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